function [T, T_deriv]=construct_chebyshev_berkson(x_interest,K,a,b,trunc_low,trunc_high,mu,sigma,fdist,alpha)

N=size(x_interest,1);
T=NaN*ones(N,K+1);

if strcmp(fdist,'Normal')==1
    for kk=1:K+1
        T(:,kk)=Tmj_Normal(x_interest,kk,0,mu,sigma,a,b);
    end
end

T_deriv=NaN*T;
for kk=1:K+1
    T_deriv(:,kk)=find_T_deriv(x_interest,0,kk,a,b,trunc_low,trunc_high,mu,sigma,fdist,alpha);
end

    %
